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ABSTRACT 

The wavelength displacement of the Diffuse Interstellar Bands at 4502, 5705, 5780, 
6284, and 7224 A with respect to the well known, narrow atomic/molecular interstellar 
lines (of Can and Nai) have been measured in the spectra of the 2 Orion Trapezium 
stars HD 37022 and HD 37020, using the HARPS-N spectrograph, fed with the 3.5 m 
Telescopio Nazionale Galileo, and the BOES spectrograph, fed with the 1.8m Korean 
telescope. The red-shift is ~25 km/s for all these DIBs. We discuss the various possible 
origins of this very peculiar wavelength shift in the light of the particular physical 
conditions in the Orion Trapezium. The above mentioned shift is seemingly absent in 
the DIBs at 6196 and 6993 A. 
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1 INTRODUCTION 

Absorption spectra of interstellar clouds contain well-known 
molecular bands of simple polar radicals, CH, CN, CH + 
which hav e been discovered and identified long ago 
dMcKellaJ Il94lh . For a long period they were believed 
to be the only possible interstellar molecules. In 1970-ties 
small homonuclear molecules (H 2 , C 2 ) were also found. 
At the same time, rotational emission features revealed 
the presence (specifically in the dense ISM, mostly in 
star-forming regions) of many complex molecules with a 
dipole momentum; the full list of these polar species (mostly 
carbon-bearing) currently counts more than 190 entries, 
clearly demonstrating that a rich chemistry, mostly car¬ 
bon-based, takes place in interstellar clouds. Some of these 
species have recently been observed in the middle and far in¬ 
frared by the space telescope Herschel. The largest currently 
known interstellar molecules are the fullerenes C 60 and C 70 
llCami et al.ll2Qlol : lBerne et a! .1120131 ). Yet unknown (uniden¬ 
tified) molecular species are believed to be the carriers of 
diffuse interstellar bands (DIBs): these make up the longest 
standing unsolved problem in astronomical spectroscopy . 
The first two DIBs were discovered in 1921 bv lHegeJ (1 1922k 
Until 1970 the list of DIBs was very short — only 9 entries. 


see e. g. http://astrochymist.org/astrochymist_ism.html 


The application of solid state detectors and more powerful 
telescopes to DIB observations led to the progressive discov¬ 
ery of many additional weak features, almost all of them in 
the range from the near-UV to the near-IR. T he current 
lists of known DIBs contain o ver 400 entries dHobbs et al.1 
l2009l : fGalazutdinov et al.ll2000k the majority of them is very 
shallow. Some of the weaker features are not unamb iguously 
detect ed and classified as DIBs in all surveys, e.g. iBondad 
J2012k Even more importantly, fine structures — reminis¬ 
cent in some cases of the unresolved rotational envelopes 
of bands of pol yatomic molecules — have been resolved in 
some DIBs (IKerr et al.ll 19981) . Nearly all conceivable forms 
of matter, from the hydrogen anion to color centres in dust 
grains, have already been proposed as DIB carriers. All of 
them have been rejected, with ve r y few still tentative excep¬ 
tions. Recently Ikrelowski et all ( 20101 1 found a very weak 
DIB which coincides with an electronic, gas-phase band of 
HC 4 H 4 ", supporting the molecular DIB origin. Another coin¬ 
cidence between a lifetime-broadened absorption spectrum 
recorded through a hydro carb on plasma and a strong DIB 
at 545 nm was reported bv lLinnartz et aD (12010k but an un¬ 
ambiguous identification of the carrier h as not yet bee n pos - 
sible. I-C 3 H 2 was initially proposed bv lMaier et al.1 1 2011 1 
but this identification would imply this species to be sur¬ 
prisingly abundant in the ISM; moreover, another broad 
DIB, initially interpreted as a second I-C 3 H 2 feature (from 
the same level) near 4883A, clearly has a different origin 
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dKrelowski et alj|201l) l. Two strong DIBs in the near-IR 
were t entatively attributed to Cg^ bv lFoing fe Ehrenfreundl 
dl994li on the basis of laboratory spectra measured in solid 
rare-g as m atrices, but this was d i spute d bv ljenniskens et al.l 
d 199711 and iGalazutdinov et al.l d2000l 'l. Despite the recent 
detec tion of the vibrational spectrum of C^ 0 (iBerne et al.l 
l2013h . its identification or rejection as DIB carrier still 
awaits g as-phase laboratory spectra, as yet unavailable. 
Similarly. llglesias-Groth et al.1 d20Qg| . l2012i ) recently claimed 
the identification of electronic transitions of gas-phase singly 
charged naphthalene, the simplest Polycyclic Aromatic Hy¬ 
drocarbon (PAH) in the spectra of Cernis 52 and, more re¬ 
cently, of HD 125241. In Cernis 52 a well-defined band is 
firmly detected, which is consistent with the strongest band 
of CioHg\ Other, weaker bands are compatible with observa¬ 
tions, but their presence is not as firmly established, mainly 
due to stellar line contamination. In the case of HD 125241, 
the authors may have been misled by the fact that the spec¬ 
trum of this star is populated with many stellar emission 
lines, possibly of circumstellar origin. One of them, centered 
at 6701.2 A and due to Siiv, is very broad. Setting the con¬ 
tinuum over the red wing of this emission can easily have 
created a broad depression as an artifact, where the authors 
identify the strongest naphthalene band. This was never con¬ 
firmed in any other sources, nor any other features of this 
molecule have been firmly identified yet on the same source, 
being at best weak features sup erimposed on stronger stel - 
lar features contaminating them. llglesias-Groth et all (|20ld l 
also claimed the identification of a band of the anthracene 
cation (CmH+q) in the spectrum of Cernis 52. However, 
this one band is obtained after subtraction of a stronger, 
very close, overlapping stellar feature, and no other bands 
of Ci 4 Hj t ~ n were detected. It was furthermore rebutted by 
IGalazutdinov et alj (1201 il l . so that this claimed identifica¬ 
tion must be regarded as tentative. No other specific PAH 
feat ures were found, despite several dedic ated efforts (see 
e. g. ISalama et al.ll201ll ; lGredel et al.ll201ll l. 

After more than 90 years of DIB research, despite 
having learned much about diffuse and translucent clouds, 
the identification of DIBs is thus still missing. We em¬ 
phasize that a vast majority of reddened stars is observed 
through several translucent clouds. As a result we ob¬ 
serve DIBs in ill-defined averages — physical parameters 
and spectra of individual clou ds may be widely different 
({Krelowski fe Westerlund Il988i l. The fact that DIBs are 
most frequently observed through a number of possibly very 
different clouds is likely to be one of the reasons for their ap¬ 
parent uniformity, washing out individual differences. This 
conversely makes all “peculiarities” very attractive, as the 
latter are most likely spectra of individual clouds, where 
the physical conditions may be more homogeneous than 
in a random sample of clouds and thus a more sensible, 
unambiguous physical interpretation of the spectra is pos- 
sible. T he examples of e xtrem e “peculiarities” shown by 
iKrelowski fe Westerlundl (1 19881 1 — namely crSco and £Oph 
— became the archetypes of translucent interstellar clouds 
of evidently different spectra. 

In this sense, objects belonging to Ori OBI asso¬ 
ciation, a region of active star formation, were discov¬ 
ered to be among the most interesting ones. The stars 
forming the famous Orion Trapezium are quite substan¬ 
tially reddened, but the clouds obscuring these objects 


are ve ry peculiar. Their extincti on curves are very un¬ 
usual dFitzpatrick fe Massal l2007tl and the features car¬ 
ried by the simple radicals typically fall below the level 
of detection. For what concerns DIBs, only two pretty 
broad ones, at 5780 and 6284 A, are strong with respect 
to extinction; the ones at 6196^ 6205, 6010, an d 5705 A 
are detectab le b ut w eak ([Krelowski fe Greenberd Il999l f . 
IKrelowski fe Greenberd dl999f l in the same paper also found 
rather surprising red-shifts of some DIBs in relation to 
the interstellar sodium doublet, the only evident atomic 
line clearly seen in the then available spect ra of the Ori 
OBI stars. This result was confirmed later by I Walker et al.l 
(l200lll . but only for HD 37061. Since the stability of the rest 
wavelengths of DIBs is one of their most well-established 
observational properties, this result is indeed extremely pe¬ 
culiar, and we therefore set out to further investigate it. 

The present paper aims to confirm the above results and 
to relate the peculiar behaviour of DIBs to that of simple 
molecular species. The point is to use high quality, high reso¬ 
lution spectra from more than one instrument. The Ori OBI 
association is the nearest region of active star formation and 
thus it is very interesting what kind of physical conditions 
may be caused by the proximity of clouds and young, hot 
stars and how these conditions may influence the formation 
of DIB carriers, DIB spectral profiles and wavelength shifts. 


2 OBSERVATIONS 

We have used the HARPS-N echelle spectrograpl0 
dCosentino et al.l l2012ll . fed with the 3.5 m Telescopio 
Nazionale Galileo. It offers a resolution of R=115,000 in 
the spectra, divided into 70 orders. This resolution allows 
to resolve Doppler components likely present in interstellar 
features, while the telescope size allows to achieve a high 
S/N ratio. Moreover, having been designed and built specif¬ 
ically for the main purpose of detecting exoplanets mea¬ 
suring radial velocity variations via accurate determinations 
of Doppler shifts, HARPS-N has a very stable and reliable 
wavelength calibration. 

In addition, to be able to rule out instrumental arti¬ 
facts, we also used data of the same targets col lected by the 
BOES echelle spectrograph ( iKim et al.ll2007l l attached to 
the 1.8m telescope of the Bohyunsan Observatory in Korea. 
This instrument can operate in three observational modes 
providing resolving powers of 30,000, 45,000 and 90,000 re¬ 
spectively. Our spectra were taken using the highest one, 
enabling us to observe relatively weak, narrow features. In 
all observational modes the spectrograph records the whole 
spectral range from ~3500 to ^10,000 A divided into 75-76 

spectral orders. _ 

All of the data were reduced using IRAF (lTodvl[l986h 
and the DECH data reduction suite, authored and made 
available^ by one of the authors (G. Galazutdinov). Data 
reduction involved the standard procedure of order location, 
bias and background subtraction, flathelding, extraction of 

1 -D spectra order by order, and wavelength calibration 
through comparison with a reference spectrum of a Th-Ar 

2 http://www.tng.iac.es/instruments/harps/ 

3 http://www.boao.re.kr/BOES/BOESppt3.files/frame.htm 

4 http://www.gazinur.com/DECH-software.html 
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Table 1. Target list with basic stellar data (from the SIMBAD database. I Wenger et alJl2000h . 


HD number 

RA 

Dec 

1 

b 

Sp. type 

V (mag) 

Eb-v 

Vrad (km/s) 

22928 

03 42 55.504 

+47 47 15.17 

150.2834 

-05.7684 

B5III C 

3.01 

0.00 

4.0 km/s 

24398 

03 54 07.922 

+31 53 01.08 

162.2891 

-16.6904 

Bllb 

2.85 

0.31 

20.6 km/s 

37020 

05 35 15.829 

-05 23 14.36 

209.0072 

-19.3853 

B0.5V 

6.73 

0.26 

28.3 km/s 

37022 

05 35 16.464 

-05 23 22.85 

209.0107 

-19.3841 

07V 

5.13 

0.31 

23.6 km/s 


lamp. The accuracy of the wavelength calibration and the 
absence of instrumental effects, which can displace the final 
wavelength solution, were verified by the cross-correlation 
of the strong telluric bands present in all studied targets 
with a reference synthetic telluric spectrum. 

The stars we observed are HD 37022 and HD 37020, 
both in the Orion Trapezium, HD 24398 as a reference “nor¬ 
mal DIBs” star, and the bright unreddened star HD 22928, 
used as a telluric line divisor. Their basic data are given in 
Table [Tj 


3 RESULTS 

HARPS-N allows very accurate and reliable determinations 
of wavelengths and this is why we used it in this project. The 
range covered by the spectrograph (from 3855 to 6910 A) 
includes two prominent doublets of Call and Nal —usu¬ 
ally very strong. The Orion Trapezium stars are also known 
to ex hibit strong nebular emission lines (see e. g. iBaldwinl 
|2000|) . apparently originating in close vicinities of the ob¬ 
served stars. This spectrograph is thus very well-suited 
to che ck w hether the re d-s hift of some DIBs reported by 
iKreiowski fc Greenberd d 19991 ) can be confirmed. 

Figure |T] shows a zoom-in of the two atomic lines of 
Nal at 5896 A and Can at 3933 A in heliocentric radial 
velocity scale, for the Orion trapezium stars HD 37020 and 
HD 37022. These stars are separated by less than 13”, and 
indeed the lines of sight towards both of them appear to tra¬ 
verse very nearly the same foreground ISM, showing almost 
identical components with the same radial velocities and 
equivalent widths. Not unusually, while both Nal and Call 
show matching components at almost the same radial ve¬ 
locities, intensity ratios are rather different, their respective 
absolute maxima are in different radial velocity components. 

Diffuse band Doppler components usually have a bet- 
ter correlation with neutr al rather than with ionized species 
dGalazutdinov et al.ll2004h , hence in the following we consid¬ 
ered the strongest Nal component as the reference velocity 
for DIBs. 

Table[2]lists the measured components, both in Nal and 
Call lines. In Call it was necessary to use two Gaussians 
around ~30 km/s to account for the slightly asymmetric 
profile of this velocity component. Only one Gaussian was 
necessary to fit the matching velocity component of Nal. The 
velocity component at ~40 km/s could only be detected in 
the Call line. The uncertainty of radial velocities was es¬ 
timated by examining the differences between values mea¬ 
sured for the same components in the two lines of the Nal 
and Call doublets respectively. 

Figure [2] unambiguously shows the red-shift of the 
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Heliocentric Radial velocity (km/s) 

Figure 1. Doppler structure of the atomic lines of Call and 
Nal lines, in radial velocity scale, for HD 37022 (top panel) and 
HD 37020 (bottom panel). In each spectrum the thick black line 
is the normalised observed flux, thin coloured lines are individual 
fitted components, the residual (difference between observation 
and fitted model) is shown below. 

5780 A DIB using the HARPS-N spectra. We used the 
HD 24398 star as a reference “normal DIBs” object. The 
latter is known as an object in which all interstellar lines 
are dominated_by a heliocentric Doppler velocity of about 
14 km/s dBondar! l2012h . with some structure barely de- 
tecta bl e only at e xtremely high resolution and S/N ratio: 
IWeltv et al.l (11996! ) resolved 8 components for Call with ra¬ 
dial velocities between -6 and 22 km/s, but with ~80% of the 
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Table 2. Heliocentric radial velocity components measured in 
the Nal line at 5896 A and Call line at 3933 A in the two Orion 
Trapezium stars HD 37022 and HD 37020. Radial velocities are 
given in km/s, and their error is less than 0.5 km/s. Equivalent 
widths are in mA. 




HD 37022 HD 37020 

rad. vel. (eq. width) 

Comp. 1 

Nal 

Call 

6.8 (48.7T1.3) 
7.0 (31.7T0.9) 

7.2 (52.6T1.8) 

8.0 (39.74T1.3) 

Comp. 2 

Nal 

Call 

19.0 (70.4T0.6) 
19.0 (25.2T0.8) 

18.7 (61.8T0.7) 
20.0 (20.8T1.1) 

Comp. 3 

Nal 

Call 

23.4 (46.0T0.6) 
22.8 (14.9T0.7) 

23.1 (49.6T0.8) 
22.9 (13.9T0.8) 

Comp. 4 

Nal 

Call 

31.9 (14.2T0.5) 
28.8 (6.7T0.3) 

31.9(13.0±0.6) 
30.6 (22.8T0.5) 

Comp. 4b 

Call 

31.6 (23.4T0.4) 

32.1 (6.4T0.3) 

Comp. 5 

Nal 

Call 

38.8 (5.9T0.6) 

38.7 (6.6T0.7) 


estimated column density being within the range from 13 to 
16 km/s. From this figure it is apparent that the maximum 
of the 5780 A DIB in HD 37022 is not found at its expected 
wavelength, but instead it is red-shifted by ~0.58 A, corre¬ 
sponding to a Doppler shift of about 30 km/s with respect 
to the radial velocity of peak absorption in Nal. While the 
Doppler structure of the Call line shows some red-shifted 
components, the red shift we measured for the 5780 A DIB 
goes clearly beyond it, where no atomic absorption com¬ 
ponent was detected. Less evidently, but measurably, the 
5780 A DIB in HD 37022 is also broader than in HD 24398, 
by about the same amount. 

In the spectrum of HD 37022 we can also detect some 
nebular emission lines, likely to arise from the Hu region 
around the same star. The Ha and two forbidden Nil lines 
seen in emission appear to have a Doppler red-shift of only a 
few km/s with respect to the strongest component of the Nal 
lines, in any case well within the velocity range spanned by 
the Doppler structure of the Call lines. The Ha appears to 
have a small Doppler shift, a few km/s, with respect to the 
Nil lines. This agrees with previous studies of th e nebular 
emission from t he Orion Trapezium llBaldwinll2000l l. Indeed, 
iBaldwinl f2000) also found much larger velocity differences 
among nebular lines due to different species in different ioni¬ 
sation states, suggesting that they arise from different parts 
of an accelerating wind. Since our observations are centered 
on the trapezium stars, fainter nebular lines are barely de¬ 
tectable in our spectra, and thus we cannot confirm this. 

The Ha line also shows a superimposed, much broader 
absorption line. Its width hinders an accurate measurement 
of its radial velocity. A similar structure, namely an emis¬ 
sion line with an underlying, much broader one in absorp¬ 
tion, can be seen in H/3. In H /3 the absorption component 
is stronger, and it can be seen to be approximately aligned 
with the radial velocity of the emission component. Upon 
checking the radial velocity of a few photospheric lines on 
our HARPS-N data, (e. g. Nm at 4003.6 A, Hei at 4026.2 A, 
Cm at 4162.9 A, Siiv at 4212.4 A, On at 4366.9 A, Om 



Figure 2. Top panel: the strongest DIB at 5780 A in the spec¬ 
trum of HD 37022, shown in radial velocity scale together with 
other interstellar and nebular lines. The zero in the velocity scale 
is chosen to coincide with the strongest component visible in the 
Di Nal line. All spectra have been normalised to the continuum, 
and shifted for clarity. Ordinate scale is arbitrary. The two dashed 
lines mark the position of the zero velocity component and the 
position of the peak of the 5780 A DIB, clearly red-shifted by 
~30 km/s. A similar plot is also shown in the bottom panel for 
the reference star HD 24398; this object shows no detectable neb¬ 
ular lines. The 5780 A DIB appears in this case (as in almost all 
cases) to be perfectly aligned (within ~2 km/s) with interstellar 
atomic features. 

at 5508.1 A, and Civ at 5812.0 A), they appear consis¬ 
tent with t he most current val ue found in the literature 
(~23.6 km/s loiivares et alfePlj l. even if there is some scat¬ 
ter between different lines, at the level of less than ~2 km/s. 
This scatter may be due to the documented variability of 
stellar lines in HD 37022, attributed to the complex inter¬ 
play ofjtsjmagnetic^eldj^st wind and some infalling mate¬ 
rial (ISim6n-Diadl2006l : IWadell2006lf . This is mostly apparent 
when comparing e.g. the Nm line at 4003.578 A and the Hei 
one at 4009.258 A the Doppler shift of the former corre¬ 
sponds to a heliocentric radial velocity of about 9.5 km/s 
while the latter of about 41 km/s. However, the vast ma¬ 
jority of the observed stellar lines is consistent with a radial 
velocity of 22-23 km/s . 

We did not include HD 37020 in Fig. [2] because it is 
almost undistinguishable from that of HD 37022, just with 
very slightly weaker DIBs and somewhat lower S/N ratio, 
hence it would not add any additional information here (see 

Fig. 0. 

Figure [3] shows a comparison of the DIBs which can 
be measured in HD 37020, HD 37022, and the reference 
star HD 24398 with HARPS-N, again in radial velocity 
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Figure 3. Five diffuse bands observed in the HARPS-N spectrum of OOri (left and middle panel) and, for reference, of £Per (right panel). 
To ease the comparison of their radial velocities, all DIBs were scaled to the same maximum intensity. All bands have been referred to 
the peak absorption in Nai lines (see Fig. n. Four of them (the relatively broad ones) show a very similar red-shift in SOri. The fifth 
one — at 6196 A — is either not red-shifted at all or the shift is below the level of detection. 


scale, using the Doppler shift of the strongest component 
of the Di Nai line in the same spectrum as reference. Four 
of the five DIBs shown, the relatively broad ones, are clearly 
red-shifted in HD 37020 and HD 37022 with respect to the 
Di Nai line marking the zero of the velocity scale. Their 
shift is compatible with a unique relative radial velocity, 
and all of them also exhibit a similar broadening with re¬ 
spect to the DIBs measured in the reference star HD 24398. 
Rather remarkably, the peak of the narrow 6196 A DIB in 
HD 37020 and HD 37022 appears to be perfectly aligned 
(within ~5 km/s) with the Doppler velocity of the atomic 
lines. The ratio of DIB intensity is also evidently different 
between HD 37022 and HD 24398: in the former the 6284 A 
DIB is stronger than the 5780 A one, while the reverse hap¬ 
pens in the latter. Morever, a weak, narrow DIB at 6287 A 
appears on the red wing of the 6284 A DIB in HD 24398, 
that is below the level of detection in HD 37022. 

To further strengthen that the observed redshift is not 
an (unlikely) instrumental artifact of HARPS-N, Fig. [4] 
shows the 5780 and 6196 A DIBs observed in HD 37022 
with BOES at slightly lower spectral resolution (90,000 
vs 115,000), which precisely match (within ~1 km/s) the 
HARPS-N data. In addition, the same Fig. [4] also shows 
two more DIBs that are out of the spectral range covered by 
HARPS-N, namely the 6993 and 7224 A DIBs. The peak of 
the former, relatively narrow, precisely matches (within ~1 
km/s) the Doppler shift of the Di Nai line, like the 6196 A 
DIB. The 7224 A DIB, conversely, a relatively broad one, 
appears again red-shifted like the 5780 A DIB. 


1.05 - 6196 \J 



0.96- 1 -■- 1 -■-l- 

-400 0 400 


Radial velocity (km/s) 

Figure 4. Four DIBs from the BOES spectra of HD 37022, in¬ 
cluding two not measurable with HARPS-N. The narrow 6993 A 
feature is evidently unshifted as the 6196A DIB, while the broad 
7224 A one shares the same red-shift of the 5780 A DIB. 


For a solid quantitative measurement of the radial ve¬ 
locity shifts visually apparent in Figs. [5] [31 and [I] we used 
the same method commonly used to precisely measure the 
radial velocities of stars, cor relat ing their observed spec¬ 
tra with a template (see e. g. lQuelodll995l . and references 
therein). We performed the correlation analysis one DIB at 
a time, using as a template the DIB profiles measured in 
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Table 3. Radial velocity displacement (in km/s) of the DIBs, ob¬ 
tained by correlating their spectra in the Orion Trapezium stars 
with the reference spectra observed towards HD 24398. Radial ve¬ 
locities for the Trapezium stars are referred to the strongest Di 
line component, assumed to mark the “interstellar rest frame”. 
For HD 24398 the “interstellar rest frame” is assumed to co¬ 
incide with the radial velocity of the molecular CH 4300 A 
line. The listed rest wavele ngths for the diffuse bands are from 
(iGalazutdinov et al.lboQOai'l . 


DIB 

HD 37022 

HD 37020 

4501.80 A 

+25T10 km/s 

+30±10 km/s 

5705.20 A 

+284=5 km/s 

+33±5 km/s 

5780.37 A 

+24±1 km/s 

+23+1 km/s 

6195.97 A 

+3±4 km/s 

+6±5 km/s 

6283.85 A 

+20±1 km/s 

+22+1 km/s 

6993.18 A 

-1+5 km/s 

n/a 

7223.91 A 

+27+1 km/s 

n/a 


HD 24398. This measurement is not coincident with veloc¬ 
ity shifts measured using DIB peaks as a reference, nor is it 
expected to be: the correlation yields velocity shifts between 
feature centroids, which would coincide with shifts measured 
on the peaks only in case of infinite signal to noise ratio and 
perfectly identical spectral shapes between spectrum and 
correlation template. Most DIBs are detectably broader in 
the Orion Trapezium with respect to HD 24398, and noise is 
clearly far from negligible. Rather remarkably, the only two 
narrow DIBs, for which the radial velocity shift thus deter¬ 
mined is most significant and least sensitive to the determi¬ 
nation method, are compatible with zero radial velocity rel¬ 
ative to the “interstellar rest frame” given by the strongest 
components in the Nal Di and CH 4300 A lines. This cor¬ 
relation procedure, as implemented both in the DECH and 
IRAF data reduction softwares, yields the radial velocities 
maximising the correlation and their standard deviations, 
via R-statistics analysis. The results are listed in Table [3j 


4 DISCUSSION AND CONCLUSIONS 

The red-shift of the DIBs at 4502, 5705, 5780, 6284, and 
7224 A observed in HD 37022 can be interpreted in two dif¬ 
ferent ways: either the peculiar physical conditions along, or 
around, HD 37022 produce intrinsically different DIB pro¬ 
files, broader and peaking at a different position, or the ob¬ 
served DIBs are produced by Doppler structure, as possibly 
hinted by all shifted DIBs sharing (or compatible with) the 
same red-sliift and broadening in velocity scale. 

To test the feasibility of this latter hypothesis, we at¬ 
tempted a conditioned least-squares fit of the red-shifted 
5780 and 6284 A DIBs in HD 37022 with a superposi¬ 
tion of two “template” DIBs, namely the ones observed in 
HD 24398. Therefore, the DIB as observed in HD 24398 was 
assumed as the “intrinsic” profile, and several copies of it 
with different scalings and Doppler shifts were used to fit the 
profile observed in HD 37022. The fit was constrained to only 
two Doppler components, one fixed at the Doppler velocity 
of the Nal lines, the other free to vary. The free parameters 
of the fit are therefore only three (for each DIB): the two in¬ 


tensities of the components (constrained to be positive) and 
the Doppler velocity of the red-shifted component. We used 
the well-known Levemberg-Marquardt method jLevenberg 
1944), as implemented in the MPFIT package (iMarkwardt 
2009|jjfor IDL. This algorithm strives to numerically reduce 
the x 2 computed between model and data, iteratively mod¬ 
ifying an initial guess for the values of the free parameters, 
until a minimum is found. It is therefore possible to fall into 
local minima, especially with noisy data, and the solution 
found will in this case depend on the initial starting point. 
When several minima exist, they will normally have differ¬ 
ent x 2 values, and some may be rejected on the sole basis of 
this. However, it may happen that more than one solution 
has an acceptable y 2 value; in this case, some of these may 
be rejected because the residual (difference between model 
and data) is clearly not white noise, i. e. with the model 
falling systematically below the data in some section of the 
spectrum and above it in another. This indeed happens with 
our fit here, with several solutions being in principle accept¬ 
able on the basis of \ 2 alone, but with only the best one (i.e. 
the one with the absolute minimum \ 2 value) producing a 
residual of white noise, with uncorrelated adjacent values. 
To make the determination of the acceptable fit more ro¬ 
bust and less subjective, we filtered both the HD 37022 and 
the template DIBs with a mild smoothing filter (Fourier 5 
points filter as implemented in the DECH data reduction 
suite). This is justified by the fact that HARPS-N spectra 
are oversampled and the DIBs being fitted are intrinsically 
rather broad. When using the filtered data, only one accept¬ 
able solution is found, which corresponds to the absolute 
minimum of \ 2 , while the unphysical ones are not even lo¬ 
cal minima any more. The resulting, very satisfactory best 
fits for the 5780 and the 6284 A DIBs are shown in Figs. [5] 
and [7] 

In both cases, the residual between the fit and the orig¬ 
inal spectrum is essentially white noise, except for the ob¬ 
vious absence of the narrow DIB at 6287 A in HD 37022, 
which is instead present in the reference HD 24398 spec¬ 
trum and thus appears as an artifact in the residual in 
Fig. [7] When slightly smoothed spectra are used for the fit, 
the solution is very robust for both the 5780 and 6284 A 
DIBs, and in both cases the additional fitted component re¬ 
sults to be red-shifted with respect to interstellar sodium 
by 34 km/s ^ v a dd < 40 km/s, the best solution being 
~35 km/s. When the same fitting procedure is applied to the 
5780 and 6284 A DIBs in HD 37020, due to lower S/N the 
solution is not unique even after smoothing, but all accept¬ 
able solutions yield nonetheless 34 km/s ^ v a dd ^ 40 km/s 
for the radial velocity of the additional component, consis¬ 
tent with what was found for HD 37022. These results can be 
compared to the radial velocity shifts listed in Table [3 How¬ 
ever, since the latter were determined by correlation with a 
template consisting of a single component, for a meaningful 
comparison one must compare them with the weighted aver¬ 
age of the radial velocities of the two components of the fit 
performed here, the weights being their relative intensities. 
This results in very nearly the same numbers of Table [3 for 
the DIBs at 5780 and 6284 A. The noise in the other DIBs 
in HD 37022 and HD 37020 is so high that the fits are not 


5 http://purl.com/net/mpfit 
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Figure 5. Upper panel: the 5780 A DIB profile in HD 37022, 
fitted with the composition of two profiles identical to that of 
HD 24398. The radial velocity difference is 35 km/s. Lower 
panel: profile difference between the “standard” reference one of 
HD 24398 and the one observed in HD 24398, both normalized 
to the same peak depth and referred to the reference frame given 
by the radial velocity of the peak absorption in the Nal atomic 
lines observed in their respective spectra. 



Figure 6. The 5780 A DIB profile in HD 37020, fitted with the 
composition of two profiles identical to that of HD 24398 (same 
as the upper panel in Fig. [5] The radial velocity difference is the 
same, namely 35 km/s. 


robust (i. e. many equivalently acceptable fit solutions are 
found) but they are compatible with the results found for 
the 5780 and 6284 A DIBs (i. e. v a dd ~ 35 km/s is always 
one of the acceptable solutions of the fit). Therefore, while 
of course we do not claim that this is the only possible in¬ 
terpretation of the data, we do find it to be possible and 
self-consistent. 

Of course, this then raises more questions: what is pro¬ 
ducing the redshifted component of the DIBs? Some resid- 



Figure 7. Same as in Fig[5] for the 6284 A DIB. The radial ve¬ 
locity difference between the two components is the same, namely 
35 km/s. Note the presence, only in HD 24398, of the narrow 6287 
DIB. 



Figure 8. The 6284 A DIB profile in HD 37020, fitted with the 
composition of two profiles identical to that of HD 24398 (same 
as the upper panel in Fig. [7] The radial velocity difference is the 
same, namely 35 km/s. 


ual, still infalling material of the young, newly-formed, very 
early-type stars? Or do they arise in the evaporating layer 
of the PDR, at the edge of the Hu region? And why is this 
additional Doppler component not detected in atomic lines? 
Upon comparing the Doppler structures of the Nal and Call 
lines in Fig. [T] and Table [2] there appears to be a tempera¬ 
ture gradient with increasing radial velocities, whereby Nal 
lines decrease while Can increases, and the highest velocity 
component is detected only in Call, but weaker. There may 
therefore be another, more red-shifted component with a 
still higher temperature, in which thus sodium and singly 
ionized calcium are absent, but this is a long stretch with 
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the available data. Observations with a substantially higher 
S/N ratio would enable us to perforin a robust fit on the 
other three red-shifted DIBs: it is clear that if 5 different 
DIBs were robustly consistent with the same Doppler struc¬ 
ture, yielding the same radial velocity as the only acceptable 
fit for the additional component of all of them, this would 
provide rather solid evidence for this interpretation. Con¬ 
versely, if each DIB were found to require a different radial 
velocity structure, this would make this interpretation just 
a pointless mathematical exercise. Furthermore, a substan¬ 
tially higher S/N ratio might also enable us to detect either 
additional Doppler components in sodium and calcium lines, 
or weak molecular features (CH, CH + , CN), which are be¬ 
low the detection limit in current data. This would be ex¬ 
tremely interesting also because there is some observational 
evidence that e. g. CH, CH + and CN features do not always 
show all the same Doppler components, with CH + features 
sometimes shifted with respect t o the other species and to 
atomic lines (see e. g. I Allenll 19941 ') . 

As to the other possible interpretation, i. e. different 
intrinsic profiles of the red-shifted DIBs in HD 37022, this 
is tempti ng in the fight o f the recent observ ations in Her- 
schef 36 (lYork et al .11201-'ll : lOka et af1l2014a l lbh . In this ex¬ 
tremely peculiar line of sight, CH and CH + were observed to 
have an anomalously high rotational temperature, and some 
DIBs were simultaneously found to be anomalously broad 
(much more than what we observe here), with a very wide 
red wing, and the peaks of some of them slightly redshifted 
(less than what we observe here). This has been tentatively 
interpreted as unresolved rotational envelopes of molecular 
bands, which in Herschel 36 appear to have an anomalously 
high rotational excitation. If this were the correct interpre¬ 
tation for HD 37022, then this would be an intermediate 
case between the “standard” excitation conditions of the 
cold neutral ISM, where “normal” DIBs arise, and those 
prevailing in Herschel 36. However, the DIB at 5797 A is 
prominent, and prominently broadened, in Herschel 36, and 
instead is undetectable in HD 37022; also, the 6196 A DIB 
is one of the broadened ones in Herschel 36, while no such 
effect is evident in our data of HD 37022 (but some lim¬ 
ited broadening might be consistent with our observations, 
hidden by the noise). Also, the 6284 A DIB is broadened 
but not redshifted in Herschel 36, while it is redshifted in 
HD 37022. To test this hypothesis, again it would be use¬ 
ful to obtain deeper observations with higher S/N, to de¬ 
tect CH, CH + and CN features, determine their rotational 
temperatures in HD 37022 and explore their relations with 
the observed profiles and positions of red-shifted DIBs. We 
remark that HD 37022 and Hersch el 36 also have_ra,ther pe¬ 
culiar extinction curves dFitzpatrick fc Massal [200^h . hint¬ 
ing that dust properties are also unusual. Actually , Her- 
sch cl 36 appears to be a tight trip le system (}Aria.s et al . 
|2010|: ISanchez-Bermudez et al.l 120141 ') . IFitzpatrick &; Massa 
( 2007 1 give for it the spectral classification of the third, most 
luminous, O 7.5 V type star of the system, which certainly 
dominates the IUE spectrum they used to derive the ex¬ 
tinction. The spectrum in the visible, in which the peculiar 
DIBs were observed, contains non-negligible contributions 
from all three unresolved stars. In any case, the extinction 
curve and DIBs are very likely to be extremely similar for 
all three stars of such a close system, unless circumstellar 
material makes a large contribution to them. 
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Figure 9. Extinction curves of HD 37020, HD 37022, 
and Herschel 36. The do tted lines are the observations by 
IFitzpatrick 8- Massal 12007 ). the shaded area represents the er¬ 
ror estimated by the authors. The contin uous lines represent the 
best fit with the [CM] 2 model ilM ill as et a] I fit) 1 .4 ) . the dashed 
lines the contribution by m acroscopic dust, th e dash-dotted lines 
the contribution by PAHs llMulas et a.l .ifooTol) . 
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Figure 10. Same as Fig. [9] but for the average Galactic extinc¬ 
tion curve and HD 23180 (o Per), another star in the same associ¬ 
ation of HD 24398 (( Per), likely to have a very similar extinction 
curve. 


Figure [9] shows the extinction curves reported by 
IFitzpatrick fc Massal (120071 ) for HD 37020, HD 37022, and 
Herschel 36, togeth er with their fits with the [CM] 2 dust 
model (IMulas et alJl2013h . Figure [TUI for comparison, shows 
the extinction curve of HD 23180 (o Per), the closest 
pr oxy of HD 24398 (C P er) available in the sample of 
iFitzpatrick fe Massal (l2007ll . and the average galactic one. 
The macroscopic dust component appears to be very sim¬ 
ilar for all the three stars in Fig. [9] all very different from 
the r eference ones i n Fig . 1101 Upon examining Tables 4 and 
5 in iMulas et al.l (120131 , one sees that in all three cases, 
upon comparison with the galactic average and HD 23180, 
the dust population is dominated by large silicatic grains, 
covered by a thick mantle of heavily processed (entirely, or 
almost entirely, aromatic) carbonaceous material. PAHs are 
very underabundant in all of them with respect to the refer¬ 
ence ones, similar to the g alactic aver a ge. In the ev olutionary 
scenario put forward by ICecchi-Pestellini et ahl d2014l) . all 
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three lines of sight appear to be composed of “old”, heavily 
processed dust. 

It is unclear whether both DIBs and the extinction 
curves are made peculiar by some common physical cause 
related to the environment, or if DIBs may be altered 
due to particular dust properties. Red shift and broad¬ 
ening in molecular electronic bands are commonly caused 
by attach m ent o f mole cules to a solid substrate (see e. g. 
iTielens fe AIlamandoIafl987l) . In an envi ronment apparentl y 
devoid of small dust particles and PAHs llMulas et al.ll2013l ) , 
in which apparently smaller nanograins and macromolecules 
are either destroyed or coalesce onto larger ones, possibly 
some DIB carriers can also partially attach to grain surfaces, 
producing what we observe in HD 37022. 

Clearly a strategy to pursue to try to find the cause of 
the red shift of some DIBs in HD 37022, besides obtaining 
better observations of HD 37022 and valid regardless of the 
interpretation, would be to find other similar cases, in order 
to understand what they have in common and thus what 
makes them different from the “standard” ones. Selecting 
lines of sight for further DIB observations choosing the ones 
with pecu lia r extincti on curves similar to that of HD 37022 
(IFitzpatrick fe Massall2007l) might be a key in this respect. 
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